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ABSTRACT In this article a new np control chart for the multiple dependent state sampling using the 
neutrosophic statistics has been introduced for the efficient monitoring of the number of defective items 
in any production process or the customer services agencies. The coefficients of the control limits of 
the proposed control chart have been determined using the neutrosophic algorithms. The efficiency of its 
quick performance has been determined by computing the neutrosophic average run lengths with respect to 
different false alarm rates under different process settings at different process shift levels. The comparison 
with the existing counterparts for the quick detecting ability of the proposed chart has also been conducted. 
The practical application of the proposed chart has been elucidated using a real-world example. 


INDEX TERMS Attribute control chart, multiple dependent state sampling, neutrosophic statistics, average 


run length, indeterminacy. 


I. INTRODUCTION 

The control chart is historic and broadly utilized the method- 
ology of the statistical process control to control the produc- 
tion of defective items or the monitoring of bad repute of any 
services providing company [1]. Quality characteristics for 
any type of goods and services can be broadly classified as an 
attribute and the variable items. Attribute quality characteris- 
tic is adopted for the production process when the produced 
items are bifurcated on the criteria as the good or defective, 
life or dead, go or not go, conforming or not conforming 
satisfy or not satisfy, etc. Firstly, the attribute quality char- 
acteristic is practically easy to use and secondly, a smaller 
sample size provides sufficient information to reach the deci- 
sion [2]. Chan et al. [3] developed a monitoring scheme for 
nonconforming items using the cumulative counts by calcu- 
lating optimal probabilities of false alarm. Epprecht, et al. [4] 
developed an adaptive attribute chart for one, two, or three 
parameters to choose an effective model. Rudisill, et al. [5] 
investigated the attribute charts using the Poisson distribution 
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for observing the smaller number of nonconforming items. 
Chiu and Kuo [6] proposed an attribute chart for multivariate 
Poisson distribution. Niaki and Abbasi [7] developed the 
multi-attribute control chart when these attributes are corre- 
lated. Lee, et al. [8] investigated the np chart using the large 
subgroup size with executed statistical constraints. 

The aforementioned chart is well-presented using the mul- 
tiple dependent state sampling (MDS) instead of the con- 
ventional single sampling scheme Hsu, et al. [9]. Due to its 
broad spectrum of the production process, it has attracted 
the attention of several researchers during the last few years. 
Conventional sampling schemes are applied only to observe 
the current status of the product whether it is defective or 
non-defective but the MDS sampling scheme keeps the his- 
tory of the product to declare the process about its produc- 
tion of defective or non-defective items. Zhou, et al. [10] 
suggested that MDS sampling improves the performance 
of the np chart because the Shewhart type charts have 
the tendency of lacking the responsiveness towards small 
changes in the production process. Reference [11] claimed 
that the efficiency of the MDS sampling plan is manifold as 
compared to a single sampling scheme as the lot sentencing 
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is based upon not only the current lot information but also 
considers the information of the already submitted lots. 
Several SPC researchers explored the MDS sampling and 
declared it an efficient scheme including Wu and Wang [11], 
Aslam, et al. [12], [13]-[20]. 

It is not very uncommon for the researchers to face the 
uncertain data in real-life situations..Smarandache [21] float 
the concept of neutrosophic statistics, which is the gener- 
alization of the traditional statistics, to analyze the uncer- 
tain, unclear, vague, and incomplete data. Smarandache [22] 
enhanced the fuzzy logic due to its limitations and sug- 
gested a neutrosophic set to deal with the doubted situa- 
tions. Khan, et al. [23] designed the dispersion chart using 
the standard deviation for the neutrosophic statistics Aslam 
and Albassam [24] presented the post hoc tested for the 
neutrosophic statistics. Albassam, et al. [25] designed a tech- 
nique to study the village population of the USA under 
the neutrosophic statistics. Zhang, et al. [26] expressed the 
general symmetry of singular neutrosophic extended group 
under the neutrosophic logic. Aslam, et al. [27] developed 
a chart for the neutrosophic statistics under the gamma 
distribution. Aslam, et al. [28] developed a new attribute 
Shewhart chart under the neutrosophic interval method. 
Abdel-Basset et al. [29] developed criteria for project selec- 
tion using a hybrid neutrosophic environment. Rich material 
is available for further reading on neutrosophic statistics for 
instance [30]-[46]. Motivated by the neutrosophic statistics 
the current article is aimed to develop an attribute np chart 
using the MDS sampling scheme. 

In this paper, we will explain an efficient monitoring 
scheme for attribute data using the neutrosophic interval 
method to determine the control chart coefficients and to 
calculate the neutrosophic average run lengths under differ- 
ent process settings. The rest of the paper is organized as 
the design of the suggested chart is explained in the sec- 
ond section. In section three the advantages of the proposed 
methodology have been discussed. The comparison of the 
proposed chart with the existing np attribute neutrosophic 
chart has been given in section four. A simulated example 
has been included and discussed in section five. 


Il. DESIGN AND METHODOLOGY OF THE 

PROPOSED CHART 

In this section, the methodology of the np chart using the 
MDS sampling under the neutrosophic statistics is elaborated. 
Let the neutrosophic probability pye [p_, py] be the neu- 
trosophic proportion of non-conforming items which fail to 
meet the specified standard criteria. Further, suppose that the 
neutrosophic sample size be denoted by nye [n_, ny]. We are 
dealing with the data consisting of uncertainty or fuzziness so 
the parameters of the attribute control chart can be determined 
using the neutrosophic statistics. Here we suppose that the 
neutrosophic probability that pyeé [pL, pu ] is the same for all 
ny€ [nL, ny] and further, suppose that the non-conforming 
numbers are denoted by Xye [Xz Xu]. According to [47] the 
above-mentioned conditions support to use the well-known 
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TABLE 1. The NARL of the proposed chart for ne[50, 50] and me[2, 4]. 









































































































































kın | [4.377,5.773] | [3.779,4.647] | [3.69,3.995] 
kon | [1.794,2.201] | [2.083,2.151] | [2.16,2.307] 
Pon | [0.202,0.211] | [0.224,0.238] | [0.21,0.221] 

c NARL NARL NARL 

0 | [200.29,204.47] | [300.37,310.48] | [380.07,376.74] 
0.01 | [188.72,190.69] | [292.55,294.69] | [354.57,344.48] 
0.02 | [176.3,176.11] | [281.12,275.7] | [327.49,311.78] 
0.05 | [137.71,132.45] | [232.28,210.03] | [245.96,220.78] 
0.08 | [102.86,95.17] | [176.76,148.33] | [176.1,150.04] 
0.10 | [83.49,75.39] | [142.99,114.93] | [138.99,114.92] 
0.20 | [28.94,23.95] | [44.55,31.44] | [42.52,31.99] 
0.30 | [11.50,9.26] | [15.46,10.64] | [15.28,11.25] 
0.40} [5.53,4.53] [6.64,4.76] [6.78,5.13] 
0.50] [3.18,2.73] [3.52,2.72] [3.66,2.94] 
0.60)  [2.13,1.93] [2.23,1.88] [2.33,2.02] 
0.70] [1.61,1.53] [1.63,1.48] [1.70,1.57] 
0.80] [1.34,1.31] [1.33,1.26] [1.38,1.33] 
0.90] [1.19,1.18] [1.18,1.14] [1.21,1.19] 
0.95| [1.14,1.13] [1.13,1.10] [1.16,1.14] 
1.00} — [1.10,1.1] [1.09,1.07] [1.11,1.1] 

TABLE 2. The NARL of the proposed chart for ne[100, 100] and me[2, 4]. 

kın | [3.562,6.14] | [5.291,6.826] | [4.934,6.559] 
Koy [2.06,2.067] [2.044,2.142] [1.939,2.303] 
Pon | [0.104,0.128] [0.114,0.117] (0.111,0.115] 

c NARL NARL NARL 

0 | [205.61,200.06] | [304.86,301.31] | [375.74,373.58] 
0.01 | [193.16,180.71] | [275.71,270.93] | [346.88,343.97] 
0.02 | [180.05,162.03] | [247.92,241.96] | [317.62,313.64] 
0.05 | [140.49,113.19] | [175.91,167.53] | [234.65,227.31] 
0.08 | [105.52,77.24] | [122.49,113.58] | [167.12,157.94] 
0.10 | [86.19,59.75] [96.04,87.52] | [132.03,122.61] 
0.20 | [31.12,18.29] [30.62,26.1] [41.63,35.92] 
0.30 | _[12.78,7.27] [11.95,9.94] | [15.58,12.99] 
0.40 | [6.27,3.74] [5.75,4.85] [7.15,5.99] 
0.50| — [3.63,2.37] [3.32,2.91] [3.94,3.41] 
0.60 | — [2.42,1.75] [2.22,2.04] [2.54,2.3] 
0.70 | [1.81,1.43] [1.68,1.61] [1.85,1.75] 
0.80)  [1.48,1.25] [1.39,1.37] [1.49,1.45] 
0.90| [1.29,1.15] [1.22,1.22] [1.28,1.28] 
0.95 | [1.22,1.11] [1.17,1.17] [1.22,1.22] 
1.00 | [1.17,1.08] [1.13,1.14] [1.17,1.17] 

















binomial probability distribution based upon the neutrosophic 
statistics which can be defined as 


n eee 
Py (xy; ny, PN) = e ) pN po, 


xn =0,1,..., nN 


() 
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TABLE 3. The NARL of the proposed chart for ne[150, 150] and me[2, 4]. 




































































kiy [3.597,4.558] [3.042,3.092] [3.49,4.8] 
T [1.927,2.084] [2.222,2.43] [2.11,2.272] 
Pon [0.108,0.144] [0.132,0.137] | [0.12,0.142] 
c NARL NARL NARL 
0 [202.79,205.27] [303.28,302.26] | [370.74,374.8] 
0.01 [185.84,187.87] [274.32,265.99] | [336.56,331.61] 
0.02 [168.5,168.73] [245.12,231.91] | [301.54,288.51] 
0.05 [119.84,112.55] [166.93,148.87] | [205.31,177.61] 
0.08 [81.91,70.24] [110.27,94.06] | [133.85,104.54] 
0.10 [63.09,50.82] [83.47,69.47] | [99.94,73.38] 
0.20 [18.44,11.87] [22.91,17.68] | [25.39,15.44] 
0.30 [6.97,4.28] [8.10,6.25] [8.56,5.14] 
0.40 [3.44,2.28] [3.75,3.07] [3.88,2.56] 
0.50 [2.13,1.58] [2.21,1.95] [2.26,1.69] 
0.60 [1.56,1.28] [1.57,1.47] [1.60,1.33] 
0.70 [1.29,1.13] [1.28,1.24] [1.30,1.16] 
0.80 [1.15,1.06] [1.14,1.12] [1.15,1.08] 
0.90 [1.08,1.02] [1.06,1.05] [1.07,1.03] 
0.95 [1.05,1.02] [1.04,1.04] [1.05,1.02] 
1.00 [1.04,1.01] [1.03,1.02] [1.03,1.01] 














Therefore, in this article, the following steps have been 
adopted to develop the np chart using neutrosophic statistics 
denoted by Nnp chart: 

Step-I Select a random sample of size nye [n,, ny] from 
any production process at each subgroup then count the defec- 
tive items Dye [DL, Dy]. 

Step-II If LCL2n < Dy < UCL2n; UCL2ne[UCL?2_, 
UCL2y] and LCL2ne [LCL2;, LCL2y ] then declare the pro- 
cess as in-control, otherwise, the process is straightforwardly 
declared as out-of-control if my subgroups are out-of-control. 

In this way, we can state that the proposed chart is the 
generalization of the traditional Shewhart np chart. The two 
control limits LCL1ne [LCL 1,, LCL 1y] is the lower control 
limit and UCLIne [UCL1,, UCL1y] is known as the upper 
control can be defined as 


CL1y = npyo + kin npo (1 — pno); 
PNo€ [Pro, puol , kine [kiL, kiu] (2) 
LCL1y = max[0, npyo — kın V npyo (1 — pwo)]; 
PNo€ [pxo, puo] , kine [kiL, kiu] (3) 
UCL2y = npno + kzy y npno (l — pno); 
Pno€ [pro, puo] , kone [k2L, kou] (4) 
LCL2y = max[0, npyno — kzn V npyo (1 — pwo)]; 


PNo€ [Pro, Puol , kon€ [k2L, k2u] (5) 
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where pyo€ [PLo, puo] is termed as the probability of the 
in-control process based upon the neutrosophic statistics and 
the term kiye [kiz, kiy] and kone [k2L, kay] are neutrosophic 
control limits coefficients. Then the probability-based upon 
the neutrosophic statistics under the in-control process at 
PN = Pnoé€ [Pro, puo] is calculated as 


P?,„Ipno = P(LCL2y < Dy < UCL2y) 
+{P(LCLIy < Dy < LCL2y) 
+P(UCL2y < Dy < UCLIn)} 

x [P(LCL2y < Dy < UCL2y)]"" (6) 


The probability of the in-control is given in Eq. (6) can be 
rewritten as 


PŇinlPNo 
= Eoo (+) pty (L — pno)" ™®™ 
E (D (+) Pury (1 = po)" 
. Seen (+) piy = pro)" | 


|UCL2y| n dy = nity 
x be ae (+) Pno (1 — pno) (7) 
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TABLE 4. Comparison of the proposed chart with Shewhart np attribute under neutrosophic statistics and traditional Shewhart np chart when n = 150 


and NARLyo = 370. 






























































Shewhart np Shewhart np attribute under 
P chart neutrosophic statistics The proposed chart 
NARL 

0 370.133 [372.021,380.802] [370.74,374.8] 
0.02 302.992 [320.879,306.17] [301.54,288.51] 
0.05 219.395 [246.195,218.741] [205.31,177.61] 
0.08 157.863 [184.106,156.677] [133.85,104.54] 
0.10 127.151 [151.004,126.101] [99.94,73.38] 
0.20 47.041 [58.374,46.673] [25.39,15.44] 
0.30 20.50 [25.634,20.216] [8.56,5.14] 
0.40 10.357 [12.816,10.085] [3.88,2.56] 
0.50 5.94 [7.19,5.694] [2.26,1.69] 
0.60 3.795 [4.464,3.584] [1.6,1.33] 
0.70 2.660 [3.03,2.481] [1.30,1.16] 
0.80 2.015 [2.224,1.866] [1.15,1.08] 
0.90 1.630 [1.747,1.507] [1.07,1.03] 
0.95 1.497 [1.584,1.386] [1.05,1.02] 
1.00 1.392 [1.455,1.292] [1.03,1.01] 











The average run length (ARL) is used to evaluate the perfor- 
mance of the developed control chart which may be defined 
as the average number of samples before the process is 
declared as out-of-control. In this article, the ARL calculated 
for the neutrosophic statistics are denoted as NARL. Thus, 
the NARL is computed as 


1 


1— Phin |PNo i 
NARLyoce [NARLyNoL, NARLyou] (8) 


NARLNo = 


The shift is considered as the natural inherent phenomena 
which may provoke any time in the running process. There- 
fore, let we suppose that the fraction of the defectives may 
exist in the process in which this fraction is shifted from 
Pn = Pno€ [pro puo] to py = pnie lpri, pui]. Then, 
the neutrosophic probability of the shifted the process when 
it is termed as in-control is written as 
PNinlPN1 = P(LCL2y < Dy < UCL2y) 
+{P(LCL1y < Dy < LCL2y) 
+P (UCL2y < Dy < UCLIy)} 


x [P(LCL2y < Dy < UCL2y)]"" (9) 
where py1 = pno + cpno; Pno€ [PLo; puol. 
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The probability Plan |pn1ı can be written as 


PNinIPN1 


| yIUCL2y _ 
ap ae |LCL2n|-+1 (> 7 Pry — pni) 
|LCL2y| hee. = 
= ae (+) Py (l — pni) 
|UCL1y| n\ oa ie 
r De ee (+) Py  ( — Pri) 


|UCL2n| n\ ag Pe my 
x pee (+) Py, C — pnt) (10) 


Then, the NARL denoted by NARLy 1 for the shifted process 
it is calculated as 


1 . 
1 — PhinlPNi 
NARLwn1€ [NARLyiz_, NARLyiuy] (11) 


NARLN = 


An R-code program was written using the above-mentioned 
methodology under different parametric settings. The control 
chart coefficients kine [kiz, kiy] and kone [k2L, koy] have 
been determined and mentioned in Tables 1-3 under different 
specified prameters. For example in Table 1 the control chart 
coefficient kıy is estimated as [4.377,5.773], [3.779,4.647] 
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TABLE 5. The simulated data for the proposed chart. 




























































































Sr# 1 2 3 4 5 6 7 8 9 10 
DN | [25,23] | [14,20] | [22,27] | [14,20] | [22,24] | [26,15] | [19,22] | [20,19] | [15,19] | [15,13] 
Sr# 11 12 13 14 15 16 17 18 19 20 
DN | [16,20] | [18,28] | [17,27] | [20,21] | [14,24] | [23,23] | [20,26] | [15,18] | [17,24] | [15,20] 
Sr# 21 22 23 24 25 26 27 28 29 30 
DN | [19,27] | [21,30] | [18,35] | [27,29] | [22,26] | [25,21] | [15,28] | [23,29] | [16,25] | [20,37] 
Sr# 31 32 33 34 35 36 37 38 39 40 
DN | [21,21] | [22,23] | [19,21] | [22,42] | [19,28] | [28,24] | [19,18] | [20,19] | [20,33] | [23,25] 
— — 
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FIGURE 1. The control chart display of the proposed chart, chart by [28] and traditional Shewhart np chart, respectively. 


and [3.693.995] and k2y is estimated as [1.794,2.201], 
[2.083,2.151] and [2.16,2.307] with pyo = [0.202,0.211], 
[0.224,0.238] and [0.21,0.221 respectively. Likewise, the 
NARL are determined for different specific NARLo = 200, 
300 and 370 and the shifted NARL are generated for several 
shift levels as 0, 0.01, 0.02, 0.05, 0.08, 0.10, 0.20, 0.30, 0.40, 
0.50, 0.60, 0.70, 0.80, 0.90, 0.95 and 1.00 under different 
samples as 50, 100 and 150 respectively. It is worth mention- 
ing here that the smaller NARL values show that the process 
is quickly announced or declared an out-of-control process. 
From Table 2, for example, the proposed chart detected the 
out-of-control process in just [3.63, 2.37], [3.32, 2.91] and 
[3.94, 3.41] samples for a shift of 0.05 when specified NARLo 
are 200, 300 and 370 respectively. 
From Table 1 through 3 it can be observed that 
1. As the shift size increases the values of NARL are 
going to decrease. 
2. As the specified NARLo increases the values NARL, 
also increases for all shift levels. 
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Ill. COMPARISON IN NARL 

In this section, the comparison of the proposed Nnp attribute 
chart for the MDS sampling using the neutrosophic statis- 
tics has been discussed. Table 4 has been generated for 
the comparison of the proposed chart with the traditional 
Shewhart np attribute chart and Shewhart np attribute chart 
based upon neutrosophic statistics proposed by [28]. From 
Table 4 it can be observed very easily that the proposed 
chart is efficient in quick detection of the out-of-control 
process. For instance, the proposed chart detects a shift of 
0.10 from 73 sample to 99® sample. The Shewhart np 
attribute under neutrosophic statistics detects the shift from 
126" sample to 151% sample. On the other hand, the tra- 
ditional Shewhart under classical statistics detects the shift 
at 127 sample. The same pattern can be observed at all 
other levels of different shifts. So it can be stated very clearly 
that the proposed chart is efficient in quickly detecting the 
out-of-control process as compared to the existing control 
charts. 
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FIGURE 2. The control chart display of the proposed chart, chart by [28] and traditional Shewhart np chart for a real data, respectively. 


TABLE 6. Neutrosophic data of juice company. 
























































Sample No. of nonconforming units D; | Sample number | No. of nonconforming 

number units 
1 [12,13] 16 [8,8] 
2 [15,15] 17 [10,10] 
3 [8,10] 18 [5,8] 
4 [10,10] 19 [13,13] 
5 [4,4] 20 [11,13] 
6 [7,7] 21 [20,20] 
7 [16,16] 22 [18,20] 
8 [9,11] 23 [24,24] 
9 [14,14] 24 [15,15] 
10 [10,10] 25 [9,12] 
11 [5,8] 26 [12,12] 
12 [6,8] 27 [7,10] 
13 [17,17] 28 [13,15] 
14 [12,15] 29 [9,9] 
15 [22,22] 30 [6,9] 














IV. COMPARISON BY SIMULATION 

The application of the proposed chart has been given 
using a simulated example in this section. We have 
used nye [150, 150], pyoe [0.1195, 0.1418], Kiyve[3.4894, 
4.7995] and K2ye[2.165, 2.2717] for shift c = 0.20 and 
the specified NARLyo = 370. The simulated values for 
these settings have been placed in Table 5. Figure | shows 
the three control charts. On observing Figure 1 it can be 
seen that the conventional chart is unable to detect the out- 
of-control situation. Shewhart np attribute under neutrosophic 
statistics by [28] shows only one observation at sample num- 
ber 32 which falls outside the upper limit. The proposed 
control chart indicates that several points are outside the 
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control limits. In addition, it can be seen that the proposed 
control chart shows out-of-control earlier than the two exist- 
ing control charts. From this simulation, it is concluded that 
the proposed chart is efficient in detecting the shift in the 
process earlier as compared to the existing control charts. 


V. APPLICATION 

In this section, we will present a real-world example for the 
illustration and the practical application of the proposed chart. 
For this purpose, we used the data of the orange juice. Accord- 
ing to [48] “Frozen orange juice concentrate is packed in 
6-oz cardboard cans. These cans are formed on a machine by 
spinning them from cardboard stock and attaching a metal 
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bottom panel. By inspection of a can, we may determine 
whether, when filled, it could possibly leak either on the 
side seam or around the bottom joint. Such a nonconforming 
can have an improper seal on either the side seam or the 
bottom panel”. The data has been shown in Table 6 for 
ready reference. For this data set the Dye [D,, Dy] values 
are plotted in Figure 2. On observing Figure 2 it can be 
seen that the proposed chart indicates that several points are 
beyond the control limits. The control chart proposed by [28] 
shows only 4 observations in the outside of the upper control 
limit. The traditional Shewhart np shows that three points 
are beyond the upper limit. From this study, it is concluded 
that the proposed control chart is the most efficient chart 
as several observations fall outside the upper control limit. 
By comparing the three charts, the proposed control chart 
clearly indicates some issues in the cans process. 


VI. CONCLUSION 

In this article, an attribute control chart using the MDS 
sampling scheme for the neutrosophic statistics was pro- 
posed. The control chart coefficients have been determined 
for different parametric settings under different shift levels. 
The proposed chart detects shift quickly as compared to the 
existing chart. The proposed chart is an efficient addition in 
the literature of the control charts for uncertain and vague 
observations. The proposed chart can further be extended for 
the multivariate case as future research. 
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